On the Uniqueness of Schwarzschild–de Sitter Spacetime

We establish a new uniqueness theorem for the three dimensional Schwarzschild–de Sitter metrics. For this, some new or improved tools are developed. These include a reverse Łojasiewicz inequality, which holds in a neighborhood of the extremal points of any smooth function. We further prove the smoothness of the set of maxima of the lapse, whenever this set contains a topological hypersurface. This leads to a new strategy for the classification of well behaved static solutions of vacuum Einstein equations with a positive cosmological constant, based on the geometry of the maximum-set of the lapse.


Introduction
A basic and fundamental question in the study of the mathematical aspects of General Relativity is the classification of the static solutions to the Einstein equations, starting from the case of vacuum solutions. The first result in this field is the celebrated Israel's Theorem [30], concerning the uniqueness of the Schwarzschild solution, among the asymptotically flat static solutions to the vacuum Einstein equations. Different proofs and refinements of Israel's result have been proposed by many authors [3,14,17,31,32,48,49,52].
A similar analysis has been performed for asymptotically hyperbolic static solutions. Wang [51] and the second author together with Herzlich [19] proved a uniqueness theorem for the Anti de Sitter solution (compare [13]), whereas Lee and Neves [41] obtained a similar result for the Kottler spacetimes with negative mass aspect (compare [26,Remark 3.4]).
In contrast with the significant amount of achievements in the case 0, very little is known when is positive, except in the locally conformally flat case [40], for perturbations of the model solutions [28] and under pinching assumptions on the curvature [4].
In order to set the stage for the problem, we recall that a (n + 1)-dimensional vacuum static spacetime is a Lorentzian manifold (X, γ ) having the form where t ∈ R, (M, g) is a n-dimensional Riemannian manifold with boundary and u : M → R is a nonnegative smooth function vanishing precisely on the boundary of M, and satisfying the vacuum Einstein equations It is clear that a static vacuum spacetime is uniquely identified by the triple (M, g, u), which will then be referred to as a static triple or static solution. The function u is usually called lapse function or static potential in the literature. The connected components of ∂ M are referred to as horizons, as they typically correspond to Killing horizons in the associated spacetime.
Since we are interested in the positive cosmological constant case, it is natural to set = n(n − 1)/2, and work under such a normalization, being clear that the general case is readily recovered through a simple rescaling. With this normalization, the most important family of n-dimensional static triples with positive cosmological constant, the so-called Birmingham-Kottler (BK) solutions (see [7,34]), can be described as where m is a real constant and ( , g ) is a (n − 1)-dimensional simply connected Riemannian manifold with Ric(g ) = (n − 2)g . (1.2) In our terminology, the Schwarzschild-de Sitter family represents a special case of BK solutions, in which the Riemannian manifold ( , g ) is taken to be a (n − 1)dimensional sphere with the canonical metric. Observe that, in order to get a static triple out of formula (1.1), it is mandatory for the lapse function u to be well defined. It is clear that this happens only if the quantity 1 − r 2 − 2mr 2−n is nonnegative for some positive values of r . This restricts the choice of the parameter m to the interval 0 m < m max , where m max is a constant only depending on n (see (5.1) for the explicit expression). All in all, we have that for any 0 < m < m max , the BK solutions (1.1) provide well defined static triples in the region {r − (m) < r < r + (m)}, where 0 < r − (m) < r + (m) are the two positive roots of 1 − r 2 − 2mr 2−n . Also, recall that the Nariai solutions (see [47]) can be defined as the limit of the BK solutions, when m → m max . An explicit expression for these is given by where again g is a Riemannian metric on satisfying (1.2).
Having the above families of model solutions in mind, the concept of virtual mass was introduced in [11,12], through a comparison algorithm (see also [6,20] for similar ideas). Building on this concept, we prove the following uniqueness result for the Schwarzschild-de Sitter black hole: For a better understanding of the above statement, we need to recall some more terminology and concepts, still given in the spirit of the comparison with the BK solutions. Depending upon the value assumed by the (normalized) surface gravity |∇u|/u max , each horizon is either of cosmological type (low surface gravity) or of black hole type (high surface gravity). When the set on which u attains its maximum separates M into two components M ± , it was shown in [11,12] how to assign a virtual mass m ± to each component; see Section 5 for more details. The virtual masses are shown there to obey the inequality m + m − , provided M + is bounded only by horizons of cosmological type (we then say that M + is an outer region). In agreement with (5.1) and by definition, for n = 3, the whole range of the virtual mass parameter m is 0 m 1/(3 √ 3). The case m = 0 leads to the de Sitter space [11,Theorem 2.3]. The case m = 1/(3 √ 3) leads to the Nariai solution. The known explicit solutions have the property that the virtual masses coincide, i.e., m + = m = m − , for some 0 < m < 1/(3 √ 3). It was shown in [12, Theorem 1.9] that, in three space-dimensions, equality of the masses and connectedness of the part of the horizon bounding M + implies that the sharp area bound (see [12,Theorem 1.4]) is saturated and thus (M, g, u) arises from the Schwarzschild-de Sitter solution. Theorem 1.1 can be seen as an improvement of this result: indeed, we are able to remove the hypothesis on the connectedness of the cosmological horizon. Perhaps more significantly, a completely new strategy of proof is used. As sketched above, the central idea in [12] was in a sense closer to the concepts developed in [3], as it consisted in establishing sharp a priori bounds relating the area of the horizons to the virtual mass of the corresponding region (Riemannian Penrose-type inequalities), characterizing the equality case in terms of the model solutions, and then showing that in some physically relevant situation (connected cosmological horizon, same mass for the outer and the inner region) the equality was actually achieved in the a priori bounds. At the heart of the new strategy is instead a detailed analysis of the locus MAX(u), established in a fairly large generality and which we believe of independent interest. Forever, before describing the most significant achievements of this analysis, let us discuss the plausibility of our assumptions, making some comments about what it is worth to expect from the level set MAX(u). On the one hand, it can be argued that containing a co-dimension 1 stratum that disconnects the underlying manifold is a highly non-generic hypothesis, even for a real analytic function like the lapse. On the other hand, it would be fair to compare our assumptions with the ones that are commonly accepted and employed to prove the classical Uniqueness Theorem for the Schwarzschild Black Hole, when = 0. In this case, it is standard to assume that the metric is asymptotically flat and that the lapse function u 0 is attaining its supremum at spatial infinity (see for example [14,Assumption 1]). Building on these two conditions, it is not difficult to realise that the blow-down limit of the large level sets of u 0 consists of a totally geodesic round sphere. Clearly, such a blow-down limit is the object that should be compared with MAX(u) in the present framework. Hence, from this point of view, it might have been natural to impose a spherical geometry on MAX(u), whereas in our Theorem 1.1 we are not even assuming that it is a smooth manifold. Analogous considerations also apply to Theorem 1.2 below, where we will suppose that MAX(u) is containing a totally umbilic component, whose induced metric obeys the Einstein equation (1.2). There is therefore a sense in which our assumptions are weaker than the ones of the classical Israel's Theorem.
Let us then describe in more detail our analysis of MAX(u). Our first main result in this context is a reverse Łojasiewicz inequality, Theorem 2.2 below, which provides an estimate on the gradient of a smooth function in terms of its increment near and away from its maximum set. This is an improved version of [12,Proposition 2.3] (see also [10,Section 1.1.5]). The result is used to control potential singularities in the function W which appears in the proof the gradient estimates (5.3) below and in turns plays a key role in the proof of Theorem 1.1. Our next main result is Corollary 3.4, which shows that the top stratum (defined in (3.1) below) of the set of maxima of an analytic function with Laplacian bounded away from zero is a smooth embedded submanifold. This result, which sheds significant new light on the structure of static solutions, is crucial for our strategy, as it allows us to invoke the uniqueness part of the Cauchy-Kovalewskaya Theorem to classify the solutions, leading to Theorem 1.2. (Geometric characterization of BK and Nariai solutions) Let (M, g, u) be a compact orientable n-dimensional, n 3, static triple with cosmological constant > 0. Suppose that the top stratum of MAX(u) is not empty and denote by one of its connected components. Then is a compact smooth hypersurface without boundary. Furthermore, if the metric induced on by g is Einstein and is totally umbilic, then the universal cover of (M, g) is isometric to either a Birmingham-Kottler or a Nariai solution.
Some comments might be in order concerning the claims in [45], seemingly rendering our Theorem 1.1 obsolete. The claims in [45] are based on calculational mistakes, which can be seen as follows: in order to justify their assertion the authors get rid of the boundary by doubling the manifold. Note that applying this to the Schwarzschild-de Sitter space metric produces a manifold with topology S 1 × S 2 . A conformal blow-up of the metric at a point using the Green function is then asserted to lead to an asymptotically flat metric with non-negative scalar curvature and zero mass. If this were correct, one could use the rigidity part of the positive energy theorem to conclude that the doubled manifold is diffeomorphic to a sphere S 3 , which is clearly not the case for S 1 × S 2 .

Summary
The rest of this paper is structured as follows: In Section 2 we will prove the reverse Łojasiewicz inequality (Theorem 2.1). In Section 3 we will analyze the regularity properties of the extremal set of an analytic function with controlled Laplacian. More precisely, in Theorem 3.1 we show how to expand a function in a neighborhood of the set of the maximum (or minimum) points. Building on this, under the hypothesis that the Hessian does not vanish, we show in Theorem 3.3 and Corollary 3.4 that the (n − 1)-dimensional part of the extremal level set is a real analytic hypersurface without boundary. It is worth remarking that the results in Sects. 2 and 3 are not exclusive of the static realm, but they hold more generally for large classes of real-valued functions. In fact, other recent applications of these same properties have appeared in [5], where they have been used to study critical metrics of the volume functional, and in [2] where the classical torsion problem is discussed. From Section 4 we start focusing exclusively on static solutions. We show that the estimates given by Theorem 3.3 allow to trigger the Cauchy-Kovalevskaya Theorem 4.1, leading to a proof of Theorem 1.2 (see Theorem 4.2). Section 5 is devoted for the most part to the statement and proof of Theorem 5.1, which is a quite general result stating that if an outer region is next to an inner region, then the virtual mass of the outer region is necessarily greater than or equal to the one of the inner region. This theorem is supplemented by a rigidity statement in the case of equality of the virtual masses. In dimension 3, it is possible to combine this rigidity statement with Theorem 1.2, leading to the proof of Theorem 1.1. Finally, in Section 6 we discuss how to exploit the Cauchy-Kovalevskaya scheme proposed in Section 4 in order to produce local static solutions.

Reverse Łojasiewicz Inequality
Let (M, g) denote a smooth n-dimensional Riemannian manifold, possibly with boundary, n 2. Given a smooth function f : M → R, we will denote by f max the maximum value of f , when achieved, and by the set of the maxima of f , when nonempty. We will assume that MAX( f ) does not meet the boundary of M, if there is one.
We start by recalling the following classical result by Łojasiewicz, concerning the behaviour of an analytic function near a critical point:  [38]) Let (M, g) be a real analytic Riemannian manifold and let f : M → R be an analytic function. Then for every point p ∈ M there exists a neighborhood U p p and real numbers c p > 0 and 1 θ p < 2 such that for every x ∈ U p it holds Let us make some comments on this result. First, we observe that the above theorem is only relevant when p is a critical point, as otherwise the proof is trivial. Another observation is that one can always set c p = 1 in (2.1), at the cost of increasing the value of θ p and restricting the neighborhood U p . Nevertheless, the inequality is usually stated as in (2.1), because one often wants to choose the optimal θ p . Let us also notice, in particular, that the above result implies the inequality in a neighborhood of any point of our manifold. The gradient estimate (2.1) has found important applications in the study of gradient flows, as it allows to control the behaviour of the flow near the critical points. The validity of the Łojasiewicz Inequality has been extended to semicontinuous subanalytic functions in [8], and a generalized version of (2.1) has been developed by Kurdyka [37] for larger classes of functions. An infinite-dimensional version of (2.1) has been proved by Simon [50], who used it to study the asymptotic behaviour of parabolic equations near critical points. A Łojasiewicz-like inequality for noncompact hypersurfaces has been discussed in [22], where it is exploited as the main technical tool to prove the uniqueness of blow-ups of the mean curvature flow. We refer the reader to [9,21,25] and the references therein for a thorough discussion of various versions of the Łojasiewicz-Simon inequality, as well as for its applications.
On the other hand, to the authors' knowledge, the opposite inequality has not been discussed yet. In this section we prove an analogous estimate from above of the gradient near the critical points. Before stating the result, let us make a preliminary observation. Suppose that we are given a Riemannian manifold (M, g) and a function f ∈ C ∞ (M). Let p ∈ M be a critical point of f . If we restrict f to a curve γ such that γ (0) = p andγ (0) = X , for some unit vector field X ∈ T p M, we have In particular, under the assumption that ∇ 2 f (X, X ) = 0 at p, we have that the left hand side of (2.2) is locally bounded. As a consequence, we immediately obtain that the inequality holds in a neighborhood of p for some constant c > 0, provided we assume that ∇ 2 f ( p)(X, X ) = 0 for every X ∈ T p M. The next theorem tells us that a slightly weaker bound is in force at the maximum (or minimum) points of f , without any assumptions on the Hessian.
Proof. The proof follows closely the one given in [12,Proposition 2.3], where the same inequality was proven for functions f satisfying (4.1). Let us start by defining the function where c > 0 is a constant that will be chosen conveniently later. We compute and taking the divergence of the above formula Fix now a connected open neighborhood of with smooth boundary ∂ . Since is compact, we can suppose that is compact as well. By definition, we have In particular, increasing the value of c if necessary, we can make F as large as desired. Since f and |∇ f | 2 are continuous and thus bounded in , and since 0 < θ < 1, it follows that, for any c big enough, we have on the whole . For such values of c, the right hand side of (2.3) is nonnegative, that is, Notice that the coefficient that multiplies w in (2.4) is negative, as f f max and 0 < θ < 1. Therefore, we can apply the Weak Maximum Principle [27, Corollary 3.2] to w in any open set where w is C 2 -that is, on any open subset of that does not intersect . For this reason, it is convenient to choose a number ε > 0 small enough so that the tubular neighborhood is contained inside , and consider the set ε = \ B ε ( ). Up to increasing the value of c if needed, we can suppose that so that w 0 on ∂ . Now we apply the Weak Maximum Principle to the function w on the open set ε , obtaining Recalling the definition of w, taking the limit as ε → 0, from the continuity of f and |∇ f |, it follows that hence we obtain w 0 on . Translating w in terms of f , we have obtained that the inequality holds in , which is the neighborhood of that we were looking for.
In particular, we have the following simple refinement, that generalizes [ Proof. From Theorem 2.2 it follows that we can choose constants α < θ < 1 and c > 0 such that and, since we have chosen θ > α, the right hand side goes to zero as we approach p. This proves the claim.

Regularity of the Extremal Level Sets
The well known Łojasiewicz Structure Theorem (established in [43], see also [36,Theorem 6.3.3]) states that the set of the critical points Crit( f ) of a real analytic function f is a (possibly disconnected) stratified analytic subvariety whose strata have dimensions between 0 and n − 1. In particular, it follows that the set where i is a finite union of i-dimensional real analytic submanifolds, for every i = 0, . . . , n−1. This means that, given a point p ∈ i , there exists a neighborhood p ∈ ⊂ M and a real analytic diffeomorphism φ : → R n such that for some i-dimensional linear space L ⊂ R n . In particular, the set n−1 is a real analytic hypersurface and is usually referred to as the top stratum of MAX( f ). In this section we show that we can get much more information about the behaviour of our function f around the maximum points that belong to the top stratum. The next theorem extends the result in [12, Proposition 2.9]. Let be a small neighborhood of p such that = ∩ MAX( f ) is contained in the top stratum and \ has two connected components + , − . We define the signed distance to as Then there is a real analytic chart (r, ϑ) = (r, ϑ 1 , . . . , ϑ n−1 ) with respect to which f admits the following expansion: Here F is a real analytic function. If we also assume that f = −ϕ( f ) for some real function ϕ, then where G is a real analytic function. Here we have denoted byφ the derivative of ϕ with respect to f , by H, h the mean curvature and second fundamental form of with respect to the normal pointing towards + , and by R, R the scalar curvatures of M and .
Remark 3.2. We have formulated Theorem 3.1 in the context of real analytic geometry because of our intended application to the classification of static solutions of Einstein equations. It is, however, clear from the proof that the conclusions of Theorem 3.1 remain true for smooth functions on smooth Riemannian manifolds, with the following modifications: first, one should assume from the outset that is a smooth hypersurface (in which case the existence of a decomposition as in (3.1) becomes irrelevant); next, the functions r , F and G and the chart (r, ϑ) are smooth but not necessarily analytic.
Proof. Let (x 1 , . . . , x n ) be a chart centered at p, with respect to which the metric g and the function f are real analytic. From the fact that p belongs to the top stratum of MAX( f ), it follows that we can choose an open neighborhood of p in M, where the signed distance r (x) is a well defined real analytic function (see for instance [35], where this result is discussed in full detail in the Euclidean space, however the proofs extend with small modifications to the Riemannian setting). More precisely, we have where φ is a real analytic function. Since r is a signed distance function, we have |∇r | = 1, which implies in particular that one of the partial derivatives of φ has to be different from zero. Without loss of generality, let us suppose ∂φ/∂x 1 = 0 in a small neighborhood of p. As a consequence, we have that the function satisfies ∂U/∂ x 1 = −∂φ/∂x 1 = 0 in . We can then apply the Real Analytic Implicit Function Theorem (see [36,Theorem 2.3.5]), from which it follows that there exists a real analytic function u : R n → R such that In other words, the change of coordinates from (r, , is real analytic. In particular f is a real analytic function also with respect to the chart (r, ϑ), where we have set Since coincides with the points where r = 0 inside , we can apply [44, Théorème 3.1] to get where A = A(r, ϑ) is a nonnegative real analytic function. For the reader convenience, we explicitly quote the invoked reference: Going on with the proof, we observe that the function f max − f achieves its minimum value 0 when r = 0, thus In particular, we can apply [44, Théorème 3.1] again and we find that where B = B(r, ϑ) is a nonnegative real analytic function. Computing the Laplacian at the points where r = 0, using the fact that the gradient of f vanishes there, we get that Applying again [44, Théorème 3.1], it follows that we can rewrite (3.3) as This concludes the first part of the proof. Notice that we could have obtained (3.4) faster by writing a priori the expansion and then using formula (3.5) below and the fact that On the other hand, this alternative approach heavily relies on analyticity, whereas we want our proof to work in the smooth setting as well, in the spirit of Remark 3.2.
We now assume that f = −ϕ( f ) and we use this to gather more information on the real analytic function F. Set ρ = {r = ρ} and observe that all ρ with ρ small enough are smooth, since (r, ϑ) = (r, ϑ 1 , . . . , ϑ n−1 ) is a real analytic chart and |∇r | = 1 = 0. In particular, of course, we have 0 = ∩ . On each ρ , the Laplacian f of f satisfies the well-known formula where n = ∂/∂r is the g-unit normal to ρ , H is the mean curvature of ρ with respect to n and f is the Laplacian of the restriction of f to ρ with respect to the metric g induced by g on ρ . Evaluating (3.5) at ρ = 0, since f = f max and |∇ f | = 0 on 0 , we immediately get Let us focus first on the quantity f . Calling g the metric induced by g on ρ and the Christoffel symbols of g , we have where the indices i, j, k vary between 1 and n − 1. On the other hand, if we assume that f = −ϕ where ϕ is a function of f , then from (3.4) we get for all i, j = 1, . . . , n − 1. From this, it easily follows that We now differentiate formula (3.5) twice with respect to r . Observing that r rr = 0 since g rr ≡ 1, it is then easy to compute that the following identities hold on the set 0 = {r = 0}: Since we also know that ∂ f /∂r = 0 and ∂ 2 f /∂r 2 = f = −ϕ( f max ) on , from the expansions above we deduce Furthermore, the classical formula for the variation of the mean curvature (see for instance [29, Lemma 7.6]) gives where in the latter equality we have used the Gauss Codazzi equation. Now that we have computed the third and fourth derivative of f , we can use this information to improve (3.4) and get the desired expansion of f . Theorem 3.1 has some interesting consequences. Let us start from the simplest one. Given a point p in the top stratum of MAX( f ) and such that f ( p) = 0, we can compute the explicit formula for the gradient of f as we approach p as Since we are assuming that |∇ 2 f | = 0 and that x is a maximum point for f , it is clear that f (x) < 0. Using the fact that eigenvalues are continuous (see for instance [33, Chapter 2, Theorem 5.1]), it follows that the eigenvalues of the Hessian of f at p are f ( p) (which is negative by hypothesis), taken with multiplicity one, and 0, taken with multiplicity n − 1. In particular, using again the continuity of the eigenvalues, restricting our neighborhood if necessary, we can suppose that the minimal eigenvalue of ∇ 2 f is simple on the whole . Since the function x → ∇ 2 f (x) is real analytic, it is known that the simple eigenvectors are real analytic in , see for instance the discussion in [33, Chapter 2, § 1] or in [39,Section 7], where a much more general statement in discussed. Therefore the vector ν extends to a real analytic unit-length vector field throughout . In particular, ν is real analytic on ∩ and the tangent space T p = ν ⊥ p is well defined. It is then easily shown (see for instance [15, p. (x 1 , . . . , x n ) centered at p such that ν = ∂/∂ x 1 in the whole . Since ∩ is smooth and ν = ∂/∂ x 1 is orthogonal to it, it follows that ∩ ⊆ {x 1 = 0}. Since f is analytic and we have shown that f max − f = 0 on an hypersurface contained in {x 1 = 0}, it is immediate to deduce that it is possible to factor out x 1 from the Taylor expansion of f max − f . It follows that f max − f = 0 on the whole {x 1 = 0}, which implies ∩ = {x 1 = 0}. In particular, is real analytic in a neighborhood of p.
To conclude that p is contained in the top stratum of MAX( f ), it remains to show that there are no other components of MAX( f ) that pass through p. In other words, it remains to show that, making smaller if necessary, we have ∩ MAX( f ) = ∩ . To this end, we first observe that, since we have already shown that is a regular analytic hypersurface around p, the proof of Theorem 3.1 can be repeated without modifications to show that formula (3.2) holds in a neighborhood of p. More precisely, there exists a real analytic chart (r, ϑ) = (r, ϑ 1 , . . . , ϑ n−1 ) such that r is the signed distance from and f has the following expansion in : Here F is real analytic. Setting C = sup |F|, we then have that in it holds that Since the points with r = 0 belong to , from the discussion above it follows that f (0, ϑ) < 0 whenever the point with coordinate (0, ϑ) belong to . If we then set the above estimate grants us that f (r, θ) < f max in \ . In particular ∩ MAX( f ) = ∩ . It follows that ∩ is contained in the top stratum , which implies in particular that p ∈ , as desired.

A Geometric Criterion for the Classification of Static Solutions
This section is devoted to the proof of Theorem 1.2. Before discussing it, let us recall briefly, following the setup carefully discussed in [12] that a static spacetime with positive cosmological constant arises as a solution of the following problem: with M compact orientable and R ≡ n(n − 1) .

(4.1)
We recall from [18,46] that a static triple (M, g, u) is necessarily real analytic, so that in particular the results discussed in Section 3 apply. As a consequence, since obviously u = −n max M (u) = 0 on MAX(u), we can apply Corollary 3.4 to deduce that the top stratum ⊆ MAX(u) is a closed real analytic hypersurface. Hence, the second fundamental form, mean curvature and induced metric are all well defined on , so that the statement of Theorem 1.2 is perfectly rigorous. In the proof, we will need the following version of the Cauchy-Kovalevskaya Theorem for systems of quasilinear PDEs; we report the complete statement for the ease of reference. We are now ready to prove Theorem 1.2, which we rewrite here for the reader's convenience. (M, g, u) be a compact orientable n-dimensional, n 3, static triple with cosmological constant > 0.

Suppose that the top stratum of MAX(u) is not empty and denote by one of its connected components. Then is a compact smooth hypersurface without boundary. Furthermore, if the metric induced on by g is Einstein and
is totally umbilic, then the universal cover of (M, g) is isometric to either a Birmingham-Kottler or a Nariai solution.
Proof. Let us start by noticing that, given any vector field X on T and called ν the normal to , from the static vacuum equations (4.1) we deduce where the latter equality follows immediately from the general expansion (3.2). We can then apply the Gauss formula to obtain On the other hand, the umbilicity of implies h = H/(n − 1) g , hence Since n 3, this implies that H is constant. Furthermore, again from expansion (3.2), together with the static equations (4.1), we get Substituting in the Gauss Codazzi equation, Let us now assume Ric = (n − 2)λg for some constant λ ∈ R. Then R = (n − 1)(n − 2) λ and we get Here we are assuming that the normal ν has been chosen so that the mean curvature H is nonnegative. As a consequence of the above formula, necessarily λ n/(n−2) and Let be a small open set with ∩ = ∅ and \ = + − , with the normal ν pointing towards + . We have already mentioned in Theorem 3.1 that the signed distance function is an analytic function in . Since clearly |∇r | ≡ 1 in , with respect to coordinates (r, ϑ) = (r, ϑ 1 , . . . , ϑ n−1 ), with ϑ 1 , . . . , ϑ n−1 coordinates on the hypersurface ∩ , we have where the coefficients g i j are functions of the coordinates (r, ϑ). Since we want to invoke the Cauchy-Kovalevskaya Theorem as stated in Theorem 4.1, it is convenient to restrict our attention to a subdomain of that can be identified with a cylinder (−ε, ε) × B n−1 . With a slight abuse of notation, we think of the coordinates (r, ϑ) as varying in (−ε, ε) × B n−1 and as a consequence, all quantities are considered as pulled back on the same domain. With respect to the coordinates (r, ϑ), observing that the normal ν coincides with ∂/∂r on , we easily see that the second fundamental form of a level set of r satisfies Formula (4.4) tells us that We have also the following initial conditions: (4.6) We would like to apply the Cauchy-Kovalevskaya Theorem 4.1 to the function by showing that f satisfies a PDE as in (4.2) coming from the equations (4.1) and the initial conditions (4.5), (4.6) on = {r = 0}. To this end, let us rewrite problem (4.1) more explicitly in terms of the derivatives of u and the metric g. We recall that, in any coordinate system, the components of the Ricci tensor satisfy where the lower order terms are polynomial functions of the components of g, the inverse of g and their first derivatives. In particular, with respect to the coordinates (r, ϑ 1 , . . . , ϑ n−1 ), the components R i j of the Ricci tensor satisfy where μ, η take the values r, 1, . . . , n − 1 whereas a, b take only the values 1, . . . , n − 1. Again the lower order terms are polynomial functions of the components of g, the inverse of g and their first derivatives. Notice that, since g is nonsingular, the inverse of g is well defined everywhere and its components g i j are analytic functions of the components g i j . Therefore, the lower order terms appearing above are analytic functions of g and the first derivatives of g. We can now rewrite problem (4.1) as i j u = lower order terms, ∂ 2 rr u + g ab ∂ 2 ab u = lower order terms, where the lower order terms appearing in (4.7) are analytic functions depending on u, g and their first derivatives. In order to apply the Cauchy-Kovalevskaya Theorem 4.1 to problem (4.7) with initial conditions (4.5), (4.6), it only remains to show that the surface is noncharacteristic for the system (4.7). In other words, we have to check that the matrix is invertible. Since u = max M (u) > 0 on , this is trivially true, hence the Cauchy-Kovalevskaya Theorem can be applied, telling us that there is a unique analytic solution to (4.7). On the other hand, the BK and Nariai solutions also solve (4.7). Moreover, it can be checked (see (5.19) and (5.13) below) that, for any λ > n/(n−2) there is a value 0 < m < m max such that the BK solution with mass m satisfies the initial conditions (4.5), (4.6) with respect to the chosen λ, whereas the Nariai solution satisfies the initial conditions (4.5), (4.6) with respect to λ = n/(n − 2) (more explicitly, on the models λ and m are related by λ = [(n − 2)m] −2/n ). It follows that our solution has to coincide with a model solution inside . This argument shows that a tubular neighborhood of is locally isometric to a tubular neighborhood of the set of the maxima of a model solution. In other words, there exists 0 < a < u max and a Riemannian covering of (M, g, u) ∩ {u > a} by a BK or Nariai solution. If a = 0, then the matrix (4.8) is invertible on the hypersurface {u = a}, so we could apply the Cauchy-Kovalevskaya Theorem again on {u = a}. It follows that the Riemannian covering can be extended up to the points where u = 0, that is, up to the horizons. Therefore, the whole (M, g, u) is isometric to a quotient of a BK or Nariai solution.

Analysis of the Interfaces and Black Hole Uniqueness Theorem
This section is mainly devoted to the proof of Theorem 5.1 below, which will be a crucial ingredient, together with Theorem 4.2, in the proof of Theorem 1.1.
We start with a quick reminder of some definitions that will be helpful in the following discussion. For more details, we refer the reader to [12]. A connected component N of M \ MAX(u) will be called region. A component S of ∂ M that belongs to the region N is called horizon of N and its surface gravity is the constant value of |∇u|/ max M (u) on S. We will distinguish different types of regions depending on the value of the surface gravity at the horizons of N . Namely, a region will be called outer, inner or cylindrical, depending on whether the maximum surface gravity of the horizons of N is smaller than, greater than or equal to √ n, respectively. Given a region N of a solution (M, g, u) of (4.1), the virtual mass μ (N , g, u) is defined as the mass of the model solution that would be responsible for the maximum of the surface gravities detected at the horizons of N , see [12,Definition 3]. The virtual mass is always a number between 0 and m max , defined as The meaning of these definitions and the interplay between virtual mass and surface gravity on the models can be appreciated in Fig. 1. It has been proven in [11,Theorem 2.3] that the virtual mass is always well defined and it is equal to zero only on the de Sitter spacetime. Finally, given two different regions A and B, if A ∩ B =: = ∅ then it follows from Corollary 3.4 that is a real analytic hypersurface with empty boundary. The hypersurface will be called interface.
We are now ready to state the following result:  μ(A, g, u)  μ(B, g, u). μ(A, g, u) = μ(B, g, u) = m then is an umbilical CMC hypersurface and the metric induced by g on has constant scalar curvature. (ii) If A is cylindrical, then B is inner or cylindrical and obviously m max = μ (A, g, u) μ (B, g, u).

(i) If A is outer, then B is necessarily inner and
then is a totally geodesic hypersurface and the metric induced by g on has constant scalar curvature.
Theorem 5.1 can be seen as a higher dimensional version of [12, Theorem 7.2 and Theorem 8.15], although the strategy of the proof proposed here is new. In the proof, we will actually compute explicitly the value of mean curvature, second fundamental form and scalar curvature of in the rigidity case (see (5.18), (5.19) and (5.20) for case (i) and (5.13), (5.14) for case (ii)). Once Theorem 5.1 is established, Theorem 1.1 follows at once (see Subsection 5.4), thanks to the geometric criterion established in Theorem 4.2.

Analytic preliminaries
One of the fundamental ingredients in the proof of Theorem 5.1 is the gradient estimate for the lapse function proven in [12,Proposition 3.3]. In Proposition 5.2 we recall that result and we outline the proof, as it represents an important application of the Reverse Łojasiewicz inequality, and more precisely of Corollary 2.3. Our aim is to show that the gradient of the lapse function |∇u| g is pointwise controlled by the corresponding quantity |∇u m | g m on the model solution of the same mass m. Here we have denoted by u m and g m the lapse function and metric of the BK solution (1.1) with mass m ∈ (0, m max ), which we recall here for the reader's convenience: In order to write a clean statement, we first need to introduce some more notation. Given m ∈ (0, m max ), we denote by the maximum value of the lapse function u m of the Schwarzschild-de Sitter solution with mass m. We then use the Implicit Function Theorem to show that, for any value m ∈ (0, m max ) and for any region N , once u is normalized so that max M (u) = u max (m), there exists a so called pseudo-radial function : N −→ R, satisfying the relationships Here r − (m) < r + (m) are defined as the radii of the two horizons of the BK solutions of mass m, whereas r 0 (m) = [(n − 2)m] 1/n represents the radius of MAX(u) in the model solution, and the sign + or − in the boundary condition ≡ r ± (m) depends on whether N is an outer or inner region, respectively. The definition of has been chosen in order to mimic the behaviour of the radial coordinate on the model solutions. Further discussions on its interpretation and well posedness can be found in [12, Subsection 2.1]. The pseudo-radial function is then exploited to prove the following fundamental gradient estimate: (M, g, u) be a solution of (4.1), and let N be an outer or inner region with virtual mass m = μ(N , g, u) < m max . Normalize u so that max M (u) = u max (m). Then it inequality holds that

Proposition 5.2. (Gradient Estimate) Let
where the function |∇u m | g m on the right hand side has to be understood as the function of one real variable that associates to some t ∈ [r − (m), r + (m)] the constant value assumed by the length of the gradient of u m on the set {|x| = t}.
Proof. The full proof of this fact has been given, with a slightly different formalism, in [12,Proposition 3.3]. Here we recall the main points. The argument leading to (5.3) is quite delicate and relies on the fact that the quantity satisfies a convenient elliptic partial differential inequality on N , namely where (g/ 2 ) and ·|· (g/ 2 ) represent the Laplacian and the scalar product of the conformally related metric −2 g. The proof of inequality (5.5) was presented in [12,Formula (3.24)] with a different formalism. The equivalence between (5.5) and [12,Formula (3.24)] is shown in the appendix. One observes that W 0 on N ∩ ∂ M by construction, since we are comparing with a model solution (M, g m , u m ) having the same mass as (N , g, u) (see [12,Lemma 2.2] for details). We now want to show that the limit lim p→MAX(u) is equal to zero. To this end, we first notice that one has an explicit formula for the gradient of the lapse function of the model solution, namely Recalling the definition of , it is easily seen that |∇u m | g m • goes to zero at the same rate of √ max M (u) − u as we approach MAX(u). Therefore, there exists a constant C > 0 such that It follows then from the Reverse Łojasiewicz Inequality (more precisely from Corollary 2.3) that W ( p) → 0, as p → MAX(u). In particular, applying the Minimum Principle on the region ε = {|W | ε} ∩ N , for every sufficiently small ε > 0, one deduces that min ε W −ε, and in turn the desired gradient estimate.
We will also need some estimate for the lapse function and the pseudo-radial function near the interface of two regions. We start by recalling the following expansion, already proven in [12, Proposition 2.6], together with a short proof for completeness. (M, g, u) be a solution to problem (4.1) and let A, B be two connected components of M \ MAX(u) with A ∩ B =: = ∅. Then, the signed distance

Proposition 5.3. Let
is an analytic function in a neighborhood of and the function u admits the expansion where H,h are the mean curvature and traceless second fundamental form of with respect to the normal ν pointing towards A.
Proof. From Theorem 3.1 we get the analyticity of the signed distance function r and the following expansion around : We also know from (4.3) that Ric(ν, ν) = 0 on . We can then apply the Gauss-Codazzi equation to obtain Substituting in the expansion for u above, we obtain (5.7).
To put our technique in a perspective, it is worth recalling that Propositions 5.2 and 5.3 were already present in [12] and have been reported for the reader's convenience. What follows from now on represents instead an original contribution of the present paper. Here, our strategy diverges from the one in [12], as the argument in [12] was based on a delicate integration by parts along the interface, whereas the following results and propositions are essentially based on the power series expansion of some relevant quantities in a neighborhood of MAX(u). The first step is indeed a detailed expansion of the pseudo-radial function about the interface between two different regions.

Proof of Theorem 5.1-(ii)
Let us start the proof of Theorem 5.1 by first addressing the cylindrical case, since it presents far less technical difficulties while containing all the main ideas. Given a solution (M, g, u) of problem (4.1), according to [12], we normalize u so that max M (u) = 1. For a cylindrical region N ⊆ M \ MAX(u), we have a gradient estimate in the same spirit of (5.3). In fact, [12,Proposition 8.2] tells us that |∇u| is bounded by the norm of the gradient of the lapse function of the Nariai solution on the corresponding level set. More explicitly, we have the following inequality: This inequality can then be employed to prove Proof. Since we have normalized u so that its maximum is 1, from (5.7) we obtain the following expansions in terms of the signed distance r : Here H is the mean curvature of with respect to the unit normal ν = ∂/∂r (which is the one pointing inside A). An easy computation now gives By definition, r is positive in A and negative in B, hence, in order for (5.11) to hold, it must be H 0, as wished.
Now we proceed to the proof of Theorem 5.1 in the cylindrical case. Since A is cylindrical, from Proposition 5.5 we get H 0, where H is the mean curvature of with respect to the normal pointing inside A. The fact that B cannot be outer follows now from Proposition 5.6 (proven in the next subsection) that tells us that if B were outer then the mean curvature of with respect to the opposite normal (the one pointing inside B) would be positive, which is a contradiction. If B is cylindrical, applying again Proposition 5.5 to both A and B, we get 0 H 0.
Therefore, H = 0, as wished. Substituting this information inside the expansions for u and |∇u| 2 , we can refine (5.12) and compute that, if H = 0 on , then From (5.11) it follows then that |h| = 0, that is, We also know from (5.3) that Ric(ν, ν) = 0 on . Substituting these pieces of information in the Gauss-Codazzi equation, we get This concludes the proof.

Proof of Theorem 5.1-(i)
We now focus on the noncylindrical case. Let (M, g, u) be a solution to problem (4.1) and consider a region with virtual mass We start by rewriting more explicitly the gradient estimate (5.3), by writing the value of |∇u m | g m • as a function of the pseudo-radial function In particular, we can rewrite (5.3) as The aim of this subsection is to compare (5.15) with the expansions discussed in Subsection 5.1 in order to deduce some consequences on the geometry of the interface . More specifically, we now prove the following analogue of Proposition 5.5: • If A is an outer region, then Since the function is clearly monotonically decreasing, necessarily we must have m A m B . Furthermore, if m A = m = m B , then formula (5.17) tells us that . (5.18) Substituting this information inside the expansions for and |∇u| 2 , we can refine (5.16) and compute that, if (5.18) holds, then From (5.15) it follows then that |h| = 0, that is, We also know from (4.3) that Ric(ν, ν) = 0 on . Substituting these pieces of information in the Gauss-Codazzi equation, we get . (5.20) This concludes the proof.

Uniqueness of the schwarzschild-de sitter black hole
We are now ready to prove Theorem 1.1, which we restate here for the reader's convenience. we deduce that is totally umbilical and that the metric g induced by g on has constant positive scalar curvature equal to 2r −2 0 (m). Since is 2-dimensional, it follows that is isometric to a sphere of radius r 0 (m) (in particular, g is Einstein).
We can then apply Theorem 1.2 to conclude that (M, g, u) is covered by a BK solution. Recall that the BK solution is a warped product of the form [r − , r + ] × S with metric dr 2 /u 2 + r 2 g S , where (S, g S ) is a simply connected 2-dimensional manifold with Ric g S = g S and u 2 = 1 − r 2 − 2 m/r . On the other hand, from the Uniformization Theorem we know that S is a 2-dimensional sphere, hence (M, g, u) must be a quotient of the Schwarzschild-de Sitter solution. The conclusion now follows easily from the fact that there are no nontrivial quotients of the sphere that preserve orientability.
Similarly, one can prove the analogous result for the Nariai solution, by employing Theorem 5.1-(ii) in place of Theorem 5.1-(i). As mentioned in the Introduction, Theorems 5.7 and 5.8 improve on the uniqueness result for the Schwarzschild-de Sitter solution proposed in [12,Theorem 7.4]. Let us briefly recall the strategy in [12]. In [12], a monotonicity argument -inspired by [3], where the case = 0 was considered -was employed to deduce the inequality between the areas of the interface and of the outer boundary. The 3-dimensional version of Theorem 5.1, and in particular the identity R = 2/r 2 0 (m), was then exploited together with the Gauss-Bonnet formula to deduce Finally, assuming connectedness of ∂ M + , the sharp and rigid Riemannian Penrose inequality was established in the same paper (see [12,Theorem 1.4]). It was then enough to put together these inequalities, to deduce that in fact they were all saturated, easily concluding the proof by rigidity.
A couple of comments are in order. First of all, this proof is markedly threedimensional as it relies heavily on the Gauss-Bonnet Theorem. Secondly, the global nature of the proof makes it so that it is necessary to add the undesirable assumption of connectedness of the outer boundary ∂ M + .
In contrast, the proof proposed in this paper is more local in nature, as it is based only on an analysis in a neighborhood of the interface. In particular, the horizons do not play any role and there is no need to require connectedness of ∂ M + . Furthermore, the proposed proof does not rely on the Gauss-Bonnet formula, but only on Theorem 5.1 and Theorem 1.2. Since both these results work in any dimension n 3, one may wonder whether it is possible to prove a higher-dimensional version of Theorem 5.7 as well.
Unfortunately the route towards such an extension seems not so easy to find, and it is even impossible from a local perspective. Indeed, we observe that the thesis of Theorem 5.1 says that is CMC umbilical with constant scalar curvature. In order to be able to invoke Theorem 1.2, one would need to further argue that this forces the metric induced on to be Einstein. In the next section we show that in general this is not the case, for dimensions n 4, if the strategy only relies on local arguments. In Proposition 6.1 we will show that for any compact, real analytic manifold and for any choice of real analytic symmetric 2-tensors g , h on , we can find a local solution to the static vacuum Einstein equations on × (−ε, ε) such that the metric induced on = × {0} is g and the second fundamental form of is h. Hence, in dimension n 4 one can clearly choose g and h in such a way that the thesis of Theorem 5.1 is satisfied (i.e., is CMC umbilical -it is enough to take h = constant · g -and the metric g induced on has constant scalar curvature), but the solution does not belong to the BK family, as long as g is not Einstein.

Local Solutions
A triple (M, g, u) will be said to satisfy the Riemannian static Einstein equations if the spacetime metric γ Lor := −u 2 dt 2 + g, ∂ t u = 0 = ∂ t g satisfies the vacuum Einstein equations, possibly with a cosmological constant. In this section we will outline how to construct such real-analytic triples (M, g, u), by mimicking the Cauchy problem in general relativity, invoking the Cauchy-Kovalevskaya theorem to justify existence. The construction is a straightforward adaptation of that of Darmois [23], a summary of the argument can be found in [16].
We note that the Cauchy-Kovalevskaya theorem, in the current setting, does not care about the signature of the metric, and that the equations for an Einstein Riemannian metric γ Riem := u 2 dt 2 + g, ∂ t u = 0 = ∂ t g, are identical to the Lorentzian ones. There is, however a slight advantage of solving the Riemannian equations instead, as then zeros of u can be smoothly mapped to axes of rotation for the metric γ Riem , avoiding the problem of equations which are potentially singular at the zeros of u. We will, however, not pursue this line of thought, as our main interest here is to solve the equations starting from a hypersurface on which u is nowhere vanishing.
Thus, we seek to construct a solution of the static Einstein equations of the form g = dr 2 + g (r ), where g (r ) is an r -dependent family of metrics on a, say compact, real analytic manifold . The initial data at r = 0 are u(0), g (0) and their first r -derivatives at r = 0, all taken to be real-analytic. The Cauchy-Kovalevskaya theorem in Gauss coordinates, as in Theorem 4.2, provides existence of an interval r ∈ (−r 0 , r 0 ) and a metric γ Lor = dr 2 −u 2 dt 2 + g (r ) =:ĝ(r ) , on {r ∈ (−r 0 , r 0 ), t ∈ R} × . The metric γ Lor will be t-independent and will satisfy the vacuum Einstein equations with a cosmological constant provided that the initial data fieldŝ g| r =0 := −u(0) 2 dt 2 + g (0) and ∂ rĝ | r =0 are chosen to be time-independent and satisfy the usual general relativistic constraint equations. Summarizing, the Cauchy-Kovalevskaya Theorem in the context of static spacetimes gives the following: Proposition 6.1. Given a compact (n − 1)-dimensional real analytic Riemannian manifold ( , g ) and real analytic functions v, w, h on with v > 0, there exists a static spacetime X = R × (−r 0 , r 0 ) × with spacetime metric γ Lor := −u 2 dt 2 + dr 2 + g (r ), where g is a collection of metrics on with g (0) = g , ∂g /∂r | = h, u | = v and ∂u/∂r | = w.
The above provides many new local solutions of the static Einstein equations. Here local refers to the fact, that the solutions might not necessarily extend to boundaries on which u vanishes.
The question then arises, which data on leads to manifolds M which are bounded by Killing horizons; equivalently, manifolds with boundary on which u vanishes. When starting from a critical level set of u, a sufficient condition for this is provided by Theorem 1.2.
Alta Matematica (INdAM) and are partially founded by the GNAMPA Project "Principi di fattorizzazione, formule di monotonia e disuguaglianze geometriche". The paper was partially completed during the authors' attendance to the program "Geometry and relativity" organized by the Erwin Schrödinger International Institute for Mathematics and Physics (ESI). PTC acknowledges the friendly hospitality and financial support of University of Trento during part of work on this paper. His research was further supported in part by by the Austrian Science Fund (FWF) under projects P23719-N16 and P29517-N27, and by the Polish National Center of Science (NCN) under grant 2016/21/B/ST1/00940.
Funding Open access funding provided by Università degli Studi di Trento within the CRUI-CARE Agreement. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/ licenses/by/4.0/.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix A Gradient estimate
In this appendix we show the equivalence between the elliptic inequality (5.5) and [12,Formula (3.24)]. To this end, we start by recovering some useful notation and formulas from [12]. First of all, we denote by ψ : [0, u max ] → R the function such that = ψ • u. The derivativeψ of ψ satisfieṡ ψ(u) = − u n−1 n − (n − 2)m and |∇u m | g m = u |ψ(u)| . (A.1) The differential inequality in [12,Formula (3.24)] is written in terms of w = u |ψ(u)| 1 − |dϕ| 2 g , whereg = g/ 2 and ϕ is the pseudo-affine function. We will not need the precise definition of ϕ but we will need the following relations involving dϕ, again proved in [12]: It is now easy to rewrite the quantity w above as w = |∇u m | g m 1 − |∇u| 2 |∇u m | 2 g m = |∇u m | g m |∇u m | 2 g m − |∇u| 2 .